AD-A058  300 


UNCLASSIFIED 


SRI  INTERNATIONAL  MENLO  PARK  CA  F/G  13/13 

CRITICAL  LOADS  FOR  REINFORCED  CONCRETE  BUNKERS. (U) 

NOV  77  A L FLORENCE  DNA001-77-C-0181 

DNA-4469F  NL 


MICROCOPY  RESOLUTION  TEST  CHAjRT 

NATIONAL  BURtAU  Of  STANDARDS  1963-.’ 


DOC  FILE  COPY;  AD  A 0 5 8 3 0 0 


4 % 


/&-B30C  2-1?- 

DNA  4469F 


CRITICAL  LOADS  FOR  REINFORCED 
CONCRETE  BUNKERS 


SRI  International 

333  Ravenswood  Avenue 

Menlo  Park,  California  94025 

November  1977 

t 

\ 

( 

Final  Report  for  Period  7 March  1977 — 31  October  1977 

CONTRACT  No.  DN A 001  -77-C-01 81 


APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  UNLIMITED. 


THIS  WORK  SPONSORED  BY  THE  DEFENSE  NUCLEAR  AGENCY 
UNDER  RDT&E  RMSS  CODE  B344077464  Y99QAXSC06222  H2590D. 


D D C 


Destroy  this  report  when  it  is  no  longer 
needed.  Do  not  return  to  sender. 


i 


' : r':  i I \ Vi  . ‘ ' !■  . AT  -£300 


O’  I 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Date  Entered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1 REPORT  NUMBER 

2.  GOVT  ACCESSIOfi-MO, 

3 RECIPIENT'S  CATALOG  NUMBER 
\ 

DNA  4469F  V 

CRITICAL  I.OADS  FOR  REINFORCED  CONCRETE  JUNKERS  , 


A.  L./ 


Florence 


- 


Final  Xep«rt.  IMInW  j 
7 Mar  IE— 31  Oct  77,  / 

'yTgwnsTWimu  imi  mtmtiftT  i/umber 

SRI  PYU-6203  ~ 

t-XniTlllfT  r.H  r.ntuT  mmaraj.. 

DNA 'jtfril - 7 7- C-/f  1 8 


9 PERFORMING  ORGANIZATION  N Alfc-M*0-A00«e55- 

SRI  International^ 

333  Ravenswood  Avenue 
Menlo  Park,  California  94025 

IT  CONTROLLING  OFFICE  NAME  AND  ADORESS  f 

Director  ( | I JJ 

Defense  Nuclear  Agency 

Washington,  D.C.  20305 

14  MONITORING  AGENCY  NAME  ft  ADDRESSf//  different  from  Controlling  Office ) 


16  DISTRIBUTION  STATEMENT  (ol  this  Report) 


10  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  UNIT  NUMBERS 

NWED  SubS£a5kyY99QAXsfc062-22 


15  SECURITY  CLAS XTo/  this 


JJJ71 


UNCLASS IFIEEN^  / 

15a.  DECLASSI  FI  CATION  rOq^NGF 
SCHEDULE  T 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  in  Block  }0,  if  different  from  Report) 


18  supplementary  notes 

This  work  sponsored  by  the  Defense  Nuclear  Agency  under  RDT&E  RMSS  Code 
B344077464  Y99QAXSC06222  H2590D. 

19.  KEY  WORDS  f Continue  on  reverse  side  if  necessary  and  identify  by  block  number) 

Nuclear  Shallow-Buried  Structures 

Blast  Pressure  Impulse 

Burst  Damage 

Bunker  Concrete 

20.  ABSTRACT  (Continue  on  reverse  side  if  necessary  and  identify  by  block  number) 

~^rAmong  the  targets  of  interest  for  nuclear  weapons  are  underground  bunkers 
of  reinforced  concrete  construction.  The  Defense  Nuclear  Agency  is  conducting 
a program  to  obtain  the  failure  loads  of  shallow-buried  bunkers  subjected  to 
air  blast  waves  acting  on  the  ground  surface.  The  blast  waves  of  interest  are 
those  from  a nearby  nuclear  weapon  of  1-kton  yield  detonated  on  or  above  the 
ground  surface. 

This  report  describes  a simple  conceptual  framework  for  intrepreting 
experimental  and  theoretical  load-damage  results.  This  framework  consists  of 

DD  . F,  ?!!”,•,  1473  EDITION  OF  I NOV  6S  IS  OBSOLETE 

^ I JAN  73  •»  • iimj’T  iccTrrrn 


/// 

/ 7 


UNCLASSIFIED _ 

SECUBITY  CLASSIFICATION  OF  THIS  PAGE  fltTirn  D«l»  Enr»r,d) 

iif*  A i u r.  r) 


• ?8  06  29 


J ( 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  P AGECHTjen  Data  Entered) 


20.  ABSTRACT  (Continued) 

a Pressure-Impulse  characterization  of  critical  loads  for  the  bunkers.  The 
procedure  for  obtaining  the  characterization  simplifies  the  prediction  of 
damage  to  box  bunkers  caused  by  ground  surface  loading.  It  requires  further 
development  for  improvement  of  the  simplifying  assumptions  and  comparisons  of 
damage  predictions  with  experimental  results. 


UNCLASSIFIED 


CONVERSION  FACTORS  FOR  U.S.  CUSTOMARY  TO 
METRIC  (SI)  UNITS  OF  MEASUREMENT 


L 


To  Convert  From To 

foot  meter  (m) 

inch  meter  (m) 

pound-force/ inch 2 (psi)  kilo  pascal  (kPa) 


3.084  x E-l 
2.540  x E-2 
6.895 


r 


TABLE  OF  CONTENTS 

LIST  OF  ILLUSTRATIONS 3 

LIST  OF  TABLES 4 

1.  INTRODUCTION  AND  SUMMARY 5 

BACKGROUND 5 

OBJECTIVE  5 

REPORT  ORGANIZATION  5 

OVERALL  RESULTS  8 

RECOMMENDATIONS  10 

2.  PRESSURE-IMPULSE  CHARACTERIZATION  OF  BUNKERS 13 

CONCEPT 13 

INTERFACE  PI  CURVE 15 

3.  SURFACE  BLAST  LOADING  17 

GROUND  SURFACE  PRESSURE  PULSE  17 

INTERACTION  PRESSURE  AND  IMPULSE 23 

THRUST  ON  ROOF  SLAB 25 

MOMENT- THRUST  YIELD  CONDITION  25 

STATIC  COLLAPSE  26 

DYNAMIC  COLLAPSE 27 

COMPARISON  WITH  WES  EXPERIMENT 29 

PRESSURE  IMPULSE  CHARACTERIZATION  31 

SUMMARY  OF  PROCEDURE 34 

4.  BURIED  CHARGE  LOADING  36 

PROCEDURE 36 

PI  CHARACTERIZATION  37 

APPENDICES 

A RESPONSE  OF  A CLAMPED  BEAM  TO  A SYMMETRICALLY  DISTRIBUTED 

RECTANGULAR  PULSE  A-l 

B REINFORCED  CONCRETE  BEAM  YIELD  CONDITION B-l 


ILLUSTRATIONS 


1.  Surface  loading  test  configuration — dynamic  load 6 

2.  Buried  charge  test  configurations  for  4-foot  and  6-foot 

ranges 7 

3.  PI  diagram  for  bunker  and  surface  loading  9 

4.  PI  predictive  method  for  buried  structure  failure  14 

5.  Rectangular  box  structure  (L/D  = 10) 18 

6.  Pressure-range  relationship  19 

7.  Impulse-range  relationship 20 

8.  Pres sure- impulse  relationship  21 

9.  Peak  pressure-pulse  duration-range  relationship  22 

10.  Pressure  pulse  20  feet  from  1-ton  surface  burst  22 

11.  Surface  loading  test  configurations — static  load 28 

12.  PI  diagram  for  bunker  and  surface  loading 33 

13.  Box  bunker  damage  from  21-pound  TNT  burst  at  2-3/4-foot 

and  4-foot  ranges 39 

14.  PI  diagram  for  bunker  and  buried  charge 40 

A1  Static  collapse  of  a rigid-plastic  clamped  beam  A-2 

A2  Dynamic  collapse  of  a rigid-plastic  clamped  beam A-5 

A3  Dynamic  response  by  initial  mechanism  A-8 

B1  Slab  geometrical  nomenclature  B-2 

B2  Stress  distributions  on  fully  plastic  cross  section  ....  B-3 

B3  Moment-thrust  yield  condition  B-ll 


3 


1. 


INTRODUCTION  AND  SUMMARY 


BACKGROUND 

An  experimental  program  is  in  progress  at  Waterways  Experimental 
Station  (WES)  to  determine  failure  loads  of  buried  reinforced  concrete 
box  bunkers.  In  the  WES  experiments,  the  bunkers  are  1/5-  or  1/10-scale 
models,  which  have  a soil  cover  of  about  half  the  bunker  height.  The 
main  type  of  loading  consists  of  high  explosive  simulation  of  ground 
surface  blast  loading  from  a nuclear  weapon  of  1-kton  yield  at  1/5-  or 
1/10-scale  (Figure  1).  The  other  type  of  loading  consists  of  a spherical 
charge  of  TNT  buried  at  bunker  midheight  and  midlength  (Figure  2). 

The  Defense  Nuclear  Agency  (DNA)  needs  a simple  conceptual  frame- 
work for  evaluating  the  experimental  results  and  for  providing  a procedure 
for  determining  critical  loads.  The  Pressure- Impulse  (PI)  characterization1 
is  being  investigated  for  this  purpose.  This  report  gives  the  status  of 
this  investigation. 

OBJECTIVE 

The  objective  of  our  work  is  to  provide  a PI  characterization  of 
critical  loads  for  the  box  bunkers.  The  procedure  for  obtaining  the 
PI  characterization  should  allow  prediction  of  the  critical  loads  for 
other  box  bunkers. 

REPORT  ORGANIZATION 

Section  2 describes  PI  characterization  and  outlines  a simple  pro- 
cedure for  approximating  the  effect  of  soil-structure  interaction. 

1G.  R.  Abrahamson  and  H.  E.  Lindberg,  "Peak  Load-Impulse  Characterization 
of  Critical  Pulse  Loads  in  Structural  Dynamics,"  Nuclear  Engineering 
and  Design,  3_7  (1),  35-46  (1976). 
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Section  3 develops  the  procedure  for  buried  box  bunkers  subjected  to 
ground  surface  blast  loading.  S rtion  4 contains  the  PI  characterization 
of  the  buried  box  bunkers  subjected  to  the  loading  from  buried  charges. 
Appendix  A presents  a dynamic  rigid-plastic  analysis  of  a clamped  beam 
to  provide  a theoretical  damage  prediction  and  an  approximate  PI 
characterization.  Appendix  B is  an  outline  of  the  derivation  of  the 
yield  condition  of  the  bunker  slab;  this  yield  condition  is  a relationship 
between  the  bending  moment  and  thrust. 

OVERALL  RESULTS 

A simple  but  very  approximate  method  has  been  devised  for  providing 
a PI  characterization  of  a buried  box  bunker  subjected  to  a ground  surface 
blast  pulse.  The  approximations  stem  mainly  from  the  representation  of 
the  soil-structure  interaction.  Because  this  phase  of  the  WES  experi- 
mental program  is  just  beginning,  only  one  static  and  one  dynamic  test  have 
been  performed;  the  dynamic  test  resulted  in  low  damage.  Thus,  a PI 
characterization  based  on  matching  the  universal  structural  PI  curve  to 
the  experimental  results  is  not  yet  possible.  The  approximate  method  in- 
corporates simple  static  and  dynamic  rigid-plastic  analyses  that  provide 
damage  predictions  and  approximations  for  the  static  collapse  pressure 
and  ideal  impulse  required  for  the  PI  characterization;  the  static  collapse 
pressure  agrees  with  the  WES  experimental  result.  The  predicted  permanent 
central  deflection  of  the  bunker  roof  in  the  WES  experiment  was  high 
(1.2  inches  compared  with  0.5-inch  central  deflection  on  a 5.6  inch-thick, 
4-foot-span  bunker  roof).  This  comparison  provides  an  overall  assessment 
of  the  assumptions  that  were  used  to  simplify  the  method.  Potential 
improvements  are  suggested  below. 

Figure  3 shows  the  PI  characterization  for  ground  surface  blast 
pulses  using  the  procedure  developed  in  Section  3.  The  results  given 
by  the  PI  characterization  and  the  WES  dynamic  experiment  show  that 


The  universal  PI  curve  is  described  in  Section  2. 
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the  critical  loads  are  highly  impulsive  and  hence  the  range  (distance 
between  surface  impact  and  the  vertical  middle  plane  of  the  bunker)  is 
very  short.  Even  for  minor  damage  this  range  is  less  than  six  times 
the  bunker  roof  or  vertical  dimension;  for  example,  at  1/10-scale  a 
2x2x8  feet  bunker  is  predicted  to  sustain  a central  roof  deflection 
of  1 inch  from  a 1-ton  nuclear  surface  burst  at  a range  of  11  feet. 

A PI  characterization  was  obtained  for  a buried  box  bunker  subjected 
to  buried  charges.  The  characterization  was  achieved  by  passing  the 
universal  structural  PI  curve  through  the  interface  PI  values  corres- 
ponding to  the  range  (distance  between  buried  charge  and  bunker  face)  for 
experimentally  observed  incipient  failure;  this  matching  determined  Iq, 

the  ideal  impulse  value  required.  The  static  collapse  pressure  P is  the 

o 

known  theoretical  value  for  a rigid-plastic  clamped  beam  subjected  to  an 
assumed  parabolic  distribution  of  interface  pressure  to  approximate  the 
effect  of  the  soil-structure  interaction. 

RECOMMENDATIONS 

These  recommendations  concern  the  PI  characterization  of  reinforced 
concrete  box  bunkers  subjected  to  ground  surface  blast  waves. 

Whenever  data  become  available  from  the  WES  dynamic  loading  experi- 
ments that  allow  at  least  an  estimate  of  the  incipient  failure  interface 
PI  values,  the  universal  structural  PI  curve  should  be  passed  through 

this  point  to  determine  the  interface  value  of  I as  described  in  Section  2. 

o 

To  determine  the  interface  static  collapse  pressure  Pq,  the  existing  results 
of  the  WES  static  loading  experiment  can  be  used.  An  an  approximation,  one 
can  assume  a parabolic  interface  pressure  distribution  that  is  zero  at  mid- 
span and  take  the  spatial  averages  to  represent  all  PI  values. 

The  theoretical  treatment  in  Section  3 of  this  report  should  be  con- 
tinued to  improve  the  assumptions,  and  an  effort  should  be  made  to  main- 
tain the  simplicity  of  the  treatment.  The  assumptions  that  should  be 
examined  concern  the  following  aspects: 
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* Blast  pulse  shapes 

* Wave  transit  times 

* Soil  wave  transmission  and  reflection 

* Soil-structure  interaction 

A ground  surface  rectangular  pulse  is  used  in  the  procedure  to 
simplify  the  analysis.  The  pressure  and  impulse  are  the  same  as  the 
blast  pulse,  which  is  more  closely  represented  by  an  exponential  decay 
following  the  instantaneous  rise  to  peak  pressure.  At  the  important 
impulsive  end  of  the  PI  curve,  the  pulse  shape  is  not  important;  hence, 
the  existing  analysis  should  be  extended  with  a rectangular  pulse  to 
treat  pressures  higher  than  5.3  times  the  static  collapse  pressure. 

Another  analytical  simplification  in  the  procedure  is  the  assumption 
that  the  roof  thrust  from  the  lateral  pressure  on  the  bunker  is  applied 
only  while  the  vertical  interface  pressure  is  acting.  For  highly  im- 
pulsive, short-duration  loads,  this  assumption  is  weak  because  of  the 
time  taken  for  the  pulse  to  pass  from  the  top  to  the  bottom  of  the  bunker; 
thus,  the  analysis  should  be  modified  to  account  for  the  wave  transit 
time.  This  step  will  introduce  the  soil  bulk  modulus  and  density  into 
the  procedure. 

Currently  we  have  not  taken  into  account  either  the  attenuation  of 
the  pulse  as  it  propagates  downward  through  the  soil  or  the  pulse  shape 
change  or  reflection  from  the  bunker  roof;  these  features  should  be 
modeled  in  a simplified  manner  for  inclusion  in  the  procedure. 

Inclusion  of  soil-structure  interaction  is  currently  modeled  as  a 
simple  redistribution  of  pressure  from  a uniform  incident  pulse  to  a 
parabolic  interface  pulse.  This  assumption  should  be  examined  for  im- 
provement. For  example,  the  WES  experiments  suggest  that  a uniformly 
distributed  impulse  followed  by  a parabolically  distributed  pressure 
pulse  may  be  more  realistic. 

We  recommend  strongly  that  a pilot  experimental  program  be  initiated, 
using  bunker  models  at  about  1/40-scale,  to  run  in  parallel  with  the  WES 
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experimental  program.  Models  of  reinforced  concrete  bunkers  at  this 
scale  are  practical  and  the  experimental  program  would  generate  data 
economically  and  rapidly,  provide  validation  for  the  PI  characterization, 
and  allow  assessment  of  the  influence  of  various  geometric  and  mechanical 
parameters;  other  bunker  shapes  can  be  fabricated  and  tested  (arches, 
for  example).  The  observations  would  be  useful  in  the  planning  of  tests 
at  larger  scales. 
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2. 


PRESSURE- IMPULSE  CHARACTERIZATION  OF  BUNKERS 


CONCEPT 

Figure  4 shows  schematically  the  construction  procedure  for  a 
pressure- impulse  (PI)  diagram.  Figures  4a  and  4b  show  the  variations  of 
ground  surface  peak  pressure  and  impulse  with  range  from  a specified 
nuclear  weapon;  Brode2*3  provides  relationships  giving  these  variations. 
Figure  4c  shows  the  "free-field"  or  surface  PI  relationship  that  results 
when  the  range  is  eliminated  from  the  plotting  coordinates;  range  becomes 
an  intrinsic  coordinate  that  is  measured  along  the  curve  and  increases 
as  the  PI  origin  is  approached.  For  buried  structures,  the  PI  loading 
is  the  interface  pressure  and  impulse.  Thus,  to  obtain  the  interface 
PI  curve  shown  in  Figure  4c  from  the  free-field  or  surface  PI  curve,  a 
procedure  is  required  that  is  based  on  the  mechanics  of  load  transfer  and 
includes  soil-structure  interaction.  Again,  range  is  an  intrinsic 
coordinate. 

The  curve  labeled  "structure"  in  Figure  4c  is  an  isodamage  curve 
for  a specified  structure  in  that  it  contains  all  combinations  of  peak 
pressure  and  impulse  for  a specified  pulse  shape  that  produce  the  same 
level  of  damage.  From  experiments  and  analyses  of  simple  structures 
subjected  to  pulse  loading,  we  have  shown1  that  the  structural  PI  curve 
has  the  general  form  shown  in  Figure  4c. 

The  interface  and  structural  PI  curves  intersect  at  a point  having 
a range  value  that  is  the  intrinsic  coordinate  on  the  interface  PI  curve. 


'H.  L.  Brode,  "Review  of  Nuclear  Weapons  Effects,”  Annual  Review  of 
Nuclear  Science,  18,  153-202  (1968). 

3 

H.  L.  Brode,  "Height  of  Burst  Effects  at  High  Overpressures,"  DASA 
Report  2056,  July  1970. 
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FIGURE  4 PI  PREDICTIVE  METHOD  FOR  BURIED  STRUCTURE  FAILURE 


If  the  isodamage  curve  is  the  PI  curve  for  failure  of  the  structure, 
the  point  of  intersection  gives  the  range  within  which  the  specified 
nuclear  weapon  will  destroy  the  structure. 

Past  work  on  PI  characterization1  has  included  an  approximate 

generalization  of  the  structural  PI  curve.  The  pressure  and  impulse 

are  usually  plotted  in  terms  of  the  ratios  P/P  and  I/I  , where  P is  the 

o o o 

static  collapse  pressure  and  I is  the  ideal  impulse  to  produce  the  same 
damage  as  that  associated  with  the  isodamage  curve.  An  approximation 
of  the  curves  of  many  simple  structures  is 

(P/P  - 1)  (1/1  - 1)  = 1 (1) 

o o 

Thus,  an  approximate  PI  curve  is  established  if  P and  I are  determined 

o o 

by  special  tests  or  approximate  analysis.  If,  during  a test  series, 

values  of  P^  and  I are  measured  that  produce  incipient  structural 

failure,  only  one  of  the  two  quantities  P and  I must  be  determined. 

o o 

If  a special  static  test  is  used  to  determine  P , Eq.  (1)  gives  for  Io> 


I 

o 


T VP  " 1 

I.  1 o 

Pn/P 
1 o 


(2) 


INTERFACE  PI  CURVE 

For  a specified  bunker,  soil,  and  burial  depth,  a series  of  experi- 
ments in  which  the  surface  PI  curve  is  used  could  provide  the  PI  character- 
ization. The  experiments  would  have  to  include  the  surface  pressure  Pq 
required  to  cause  static  failure  of  the  bunker  and  a pair  of  values  P , 

1^  for  the  surface  pulse  that  causes  dynamic  failure  of  the  bunker; 
the  value  of  Iq  could  then  be  calculated  by  (2).  However,  this  procedure 
assumes  that  the  structural  PI  curve  given  by  (1)  applies  to  buried 
structures  and  it  provides  neither  a method  of  generalizing  the  result 
to  other  bunker  designs  nor  an  understanding  of  the  mechanisms  involved. 
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Thus,  a procedure  is  required  to  generate  an  interface  PI  curve  that 
accounts  for  the  mechanisms  of  load  transfer  from  the  surface  to  the 
bunker . 

The  load  transfer  process  consists  of  transmission  of  the  surface 
pulse  downward  through  the  soil,  reflection  and  diffraction  from  the 
bunker,  and  the  resulting  soil-structure  interaction.  This  complex 
process  transforms  the  fairly  uniform  surface  pulse  into  an  interaction 
pulse  having  a different  peak  pressure,  shape,  and  distribution.  If  a 
general  procedure  is  to  be  formulated  for  determining  the  interface  PI 
curve  from  the  free-field  PI  curve,  we  must  introduce  simplifying 
assumptions  to  reduce  the  complexity  of  the  load  transfer  mechanics. 
These  assumptions  must  lead  to  results  that  agree  approximately  with 
experimental  results.  Also,  the  procedure  requires  a choice  of  defini- 
tion of  the  interface  pressures  and  impulses  because  of  the  spatial 
variations.  The  development  and  application  of  the  procedure  is  the 
subject  of  Section  3. 
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3. 


SURFACE  BLAST  LOADING 


We  develop  here  the  procedure  to  obtain  a Pressure-Impulse  (PI) 
characterization  of  a shallow-buried  box  bunker  subjected  to  surface 
blast  loading.  The  procedure  is  summarized  at  the  end  of  this  section. 
Calculations  are  based  on  the  design  shown  in  Figure  5 which,  for  the 
purpose  of  illustrating  the  procedure,  is  assumed  to  be  1/10  scale. 
Calculations  are  also  based  on  a 1 kton  surface  burst  reduced  to  1/10  scale 
to  correspond  to  the  assumed  scale  of  the  box  bunker  design. 

GROUND  SURFACE  PRESSURE  PULSE 

ry 

Approximate  formulas^  relating  the  peak  pressure,  impulse,  range, 
and  yield  (P,  I,  R,  W)  are  given  in  psi,  psi-sec,  kft,  and  Mton  units, 
respectively,  are 

1 1/23/? 

P = 3300  W/R  + 192  W-  ' /R ' 

1/2  1/3  1/2 

I = 1.83  P W 7 [1  + 0.00385P  ] 

Figures  6,  7,  and  8 show  the  pressure-range,  impulse-range,  and 
pressure- impulse  curves  that  result  from  formulas  (3)  and  (4)  for  a 
surface  burst  of  1-ton  yield,  which  is  the  1/10-scale  yield  for  a 1-KT 
full-scale  yield.  The  results  in  Figures  6,  7,  and  8 are  in  approximate 
agreement  with  the  corresponding  results  from  reference  3.  A plot  of 
the  pulse  duration3  (period  of  positive  overpressure)  is  shown  in  Figure  9. 
An  example  of  the  pulse  shape3  is  shown  in  Figure  10. 

| Analysis  of  the  dynamic  response  of  rigid-plastic  structures  to 

blast  loading  can  be  greatly  simplified  by  replacing  the  blast  pulse 
by  a rectangular  pulse  with  the  same  peak  pressure  and  impulse;  by  this 
method,  however,  damage  is  overpredicted  for  pulses  with  durations  com- 
parable to  the  fundamental  structural  period;  for  rigid-plastic  structures, 
this  fundamental  structural  period  may  be  taken  as  the  duration  of  motion. 
However,  calculations  presented  below  indicate  that  durations  of  pulses 
causing  failure  are  much  shorter  than  durations  of  motion  and  are  there- 
fore impulsive. 
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INTERACTION  PRESSURE  AND  IMPULSE 


As  trial  approximations,  the  effects  of  the  soil-structure  inter- 
actions are  taken  as:  (I)  a redistribution  of  pressure  from  a uniform 
distribution  on  the  ground  surface  to  a parabolic  distribution  on  the 
bunker  roof,  having  the  minimum  pressure  at  the  center  and  the  maxima 
at  the  edges,  and  (2)  a preservation  of  the  pulse  duration  and  pulse 
shape.  The  average  pressure  of  the  parabolic  pressure  distribution  is 
equated  with  surface  pressure,  as  is  the  pulse  duration.  The  influence 
of  the  blast  wave  transit  time  is  neglected. 

From  the  infinite  possibilities  of  parabolic  distributions  of 
pressure,  we  have  chosen  the  one  producing  no  interaction  pressure  at 
the  roof  center  because  it  maximizes  the  load  required  to  cause  a 
specified  damage. 

For  the  box  design  shown  in  Figure  5,  which  has  internal  dimensions 
of  24  x 24  inches  end  a wall  thickness  h = 2.8  inches,  the  parabolic 
pressure  distribution  is 

P ■ 12(f)  «> 

where  i = 12  inches  is  regarded  as  the  half-span  of  a clamped  beam 
modeling  the  roof.  The  average  pressure  on  the  entire  top  surface  of 
length  2SL  + 2h  is 

i+h 

p = i+h  f p(x)dx  = 3 q2  (nr) 

o 

and,  by  assumption,  equals  the  corresponding  uniform  ground  surface 
pressure.  The  average  pressure  on  the  clamped  beam  of  length  2£  is 

l 

Pd  = i / P(x>dx  = 3 q2  (6) 

o 
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Hence,  the  interaction  pressure  required  in  the  PI  characterization  is 


which,  for  the  bunker  of  Figure  5,  gives 


Preservation  of  the  pulse  shape  and  duration  gives 


(7) 


(8) 


(9) 


As  a numerical  example,  Table  1 gives  the  surface  pressure  and  impulse 
at  a range  of  20  feet  as  P = 479  psi,  and  I = 0.414  psi-s  so  that,  by 
Eqs.  (8)  and  (9),  the  interaction  pressure  and  impulse  are  P^  = 319  psi 
and  1^  = 0.276  psi-s.  The  pressure  over  the  clamped  beam  supports  at 
x = ±1  is  q^  = 3p^  = 957  psi;  over  the  edge  of  the  bunker  at  x = ±1  (£  + h), 
the  pressure  is  1456  psi. 


TABLE  1 PULSE  DATA  FOR  1-TON  SURFACE  BURST 


Equivalent 


Range 

(ft) 

Peak 

Pressure 

(psi) 

Impulse 

(psi-s) 

Pulse 

Dura  tion 
(ms) 

Rectangular 
Pulse  Duration 
(ms) 

15 

1064 

0.623 

10.11 

0.586 

20 

479 

0.414 

9.11 

0.864 

25 

266 

0.297 

8.45 

1.117 

30 

168 

0.226 

8.21 

1.345 

35 

114 

0.182 

8.32 

1.597 

40 

83 

0.153 

8.71 

1.843 
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THRUST  ON  ROOF  SLAB 


In  the  soil  away  from  the  side  of  the  bunker,  the  horizontal  and 


vertical  compressive  stresses  and  become  the  principal  stresses. 


If  these  stresses  cause  distortional  yielding  of  the  soil  according  to 
the  Mohr-Coulomb  criterion,  they  satisfy  the  relationship 


o = 

X 


1 - sin4>  2c  cosj 

1 + sin<{  y 1 + sin<J> 


(10) 


where  <j>  and  c are  the  friction  angle  and  cohesive  strength  of  the  soil, 
respectively.  Away  from  the  bunker  side,  we  assume  that  the  place  wave 
descending  into  the  soil  is  replaced  by  a pressure  pulse  acting  in- 


stantaneously along  the  full  height  of  the  bunker  so  that  = P(t), 


the  surface  pressure  pulse.  The  other  stress  component  is  given  by 


Eq.  (10).  This  uniformly  distributed  lateral  pressure  is  redistributed 
on  the  side  of  the  bunker  because  cf  the  structural  response,  but  the 


average  pressure  is  taken  to  be  and  is  assumed  to  act  as  long  as  the 


roof  is  being  loaded  from  above.  The  thrust  produced  by  this  quasi-static 


model  is  therefore  N = (£  + h)a  . The  consequences  of  neglecting  wave 


propagation  are  not  considered  in  this  analysis. 

For  soil  with  <j>  = 35°  and  c = 10.7  psi,  Eq.  (10)  becomes 


a = 0.271a  - 11 

x y 


(11) 


For  a surface  pressure  of  = P = 479  psi,  we  have  = 119  psi;  hence, 

2.8  in.)  receives  a roof  thrust 


the  bunker  of  Figure  5 (£  = 12  in.,  h 
of  N = 3758  pounds. 


MOMENT-THRUST  YIELD  CONDITION 


In  the  analysis  of  Appendix  B,  Figure  B3  shows  the  moment-thrus t 
yield  condition  for  the  bunker  roof  slab.  The  data  used  are  as  follows: 


Slab  depth 

Reinforcement  spacing 
Upper  steel  depth 


h = 2.8  in. 
b = 2.0  in. 
d1  = 0.4  in. 
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Lower  steel  depth 
Diameter  of  steel  bars 
Area  of  steel  bars 
Steel  modulus 
Concrete  modulus 
Modular  ratio 
Steel  yield  strength 
Concrete  crush  strength 
Strength  ratio 
Steel  area  ratio 

The  fully  plastic  moment  and  th 

be 


d = 2.4  in. 

d = 0.25  in. 
s 

A = trd2/4  = 0.049  in.’' 
s s 

E = 28  x 106  lb/in.2 
s 

E = 3.5  x 106  lb/in. 2 
c 

n = E /E  =8 
s c 

a = 72,000  psi 

y 

a = 6,000  psi 
u 

S = a /a  =12 

y u 

a = A /bh  = 0.00877 
s 


t per  unit  slab  width  were  found  to 


M = 4102  lb. in. /in.  N = 20,042  lb. /in. 
o o 

An  approximation  to  that  portion  of  the  yield  condition  for 
0 < N/N_ < 0.1  is 

0<  < 0.1  (12) 

N 

o 


M_  = 
M 


1 + 


. N 
’iT 


STATIC  COLLAPSE 


The  parabolic  pressure  distribution  (see  Eq.  5)  applied  statically 
to  a clamped  rigid-plastic  beam  causes  collapse  in  a three-hinge 
mechanism  when  the  average  pressure  becomes 


P 

s 


(13) 


If  we  neglect  thrust  initially  and  use  the  fully  plastic  moment  Mq  in 
Eq.  (13),  the  static  collapse  pressure  is  Pg  = 228  psi.  Equation  (7) 
relates  the  average  interaction  pressure  to  the  surface  pressure  P,  so 
that 

„ l>.+h\  P = “ P *342  psi 

p-nr)  s 2 ” 

as  a first  approximation. 
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Inserting  this  value  for  a in  Eq.  (11)  gives  the  average  pressure 
on  the  side  of  the  bunker  as  = 83  psl,  resulting  in  a thrust  on  the 
roof  slab  of  N = 1,228  lb. /in.  By  Eq.  (12),  this  thrust  (N/Nq  = 0.061) 
increases  the  resistive  moment  to  M = 1.31  M^,  the  average  interaction 
pressure  to  Pg  = 298  psi,  and  the  surface  pressure  to  P = 446  psi,  as 
a second  approximation.  Three  further  iterations  give  N/N  = 0.10, 

M/M^  = 1.51,  Pg  = 344  psi,  and  P = 516  psi.  This  surface  pressure  com- 
pares reasonably  well  with  the  results  of  the  WES  static  test  conducted 
with  the  configuration  of  Figure  11. 

DYNAMIC  COLLAPSE 

Appendix  A presents  an  analysis  of  the  response  of  a clamped 
rigid-plastic  beam  to  a rectangular  pulse  with  a parabolic  distribution 
symmetric  about  the  center;  the  pressure  is  zero  at  the  beam  center. 

The  results  provide  a theoretical  prediction  of  the  central  deflection 
and  duration  of  motion  for  comparison  with  the  results  of  the  WES 
dynamic  experiment.  The  analysis  also  provides  guidance  for  the 
selection  of  an  impulse  to  play  the  role  of  I in  a PI  characterization; 
the  analysis  is  not  complete,  but  it  treats  pressures  up  to  five  times 
the  static  collapse  pressure. 

The  final  central  deflection  and  motion  duration  formulas  obtained 

are 

Wf  - ^ + A X"sPi2{2  - (l  - pi)2jJ 

~ -1+f  )»soj{3  - (l  - Pj)2} 

where 

(i  - 5^(3  + 6PX  + 5p2)  - ^ 


(15) 


(16) 
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The  range  of  A is  1 < A < 5.3  for  the  structure  of  Figure  5.  The 
symbols  denote 

W^.:  central  deflection 

p^:  average  pressure  on  beam 

t^:  rectangular  pulse  duration 

p:  beam  mass  per  unit  length 

P^:  dimensionless  initial  location  of  plastic  hinge 

(p^  = x^/£  , where  H is  the  half-span) 

A:  pressure  ratio  (A  = 'p^/p'  , where  'p  is  the  static  collapse 

pressure) 

m : static  moment  ratio  (m  = M /M  , where  M and  M are  the 
s s s o so 

yield  moments  at  static  collapse  with  and  without  thrust) 

m,:  dynamic  moment  ratio  (m , = M,/M  , where  M,  and  M are  the 

d ddo  d o 

yield  moments  during  dynamic  collapse  with  and  without  thrust), 
t^  = beam  motion  duration. 

COMPARISON  WITH  WES  EXPERIMENT 

The  WES  dynamic  test,  Essex  Box  HEST  1 (9  August  1977)  was  per- 
formed on  a reinforced  concrete  box  having  the  internal  dimensions  of 
4 x 4 x 16  feet.  For  our  comparison,  the  WES  data  are  scaled  to  apply 
to  a bunker  with  internal  dimensions  of  2 x 2 x 8 feet.  The  box  design 
is  now  given  by  that  shown  in  Figure  5,  which  we  again  assume  to  be 
1/10  scale. 

The  deflection  predicted  is  obtained  as  follows: 

* Measured  peak  surface  pressure  (WES  gage  BP3)  , P = 1800  psi 

* Measured  surface  impulse,  I = y x 3.6  = 1.8  psi-s 

* Duration  of  equivalent  rectangular  pulse,  t^  = I/P  = 1 ms 

* From  Eqs.  (8)  and  (9),  the  average  interface  pressure  and 
impulse  are 

p’j  = 1200  psi  Td  = 1.2  psi-s 

* Average  interface  static  collapse  pressure,  ps  = 344  psi 
(This  pressure  corresponds  to  a static  surface  pressure  of 
P = 516  psi) 

* Pressure  ratio  A = p’./p’  = 3.49,  so  1 < A < 5.3. 

d s 
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* Initial  location  of  plastic  hinges,  from  (16),  is  P = 0.61 

* Central  deflection,  from  (14),  is  W^.  = 2.06  in. 

* Duration  of  motion,  from  (15),  is  t^.  = 4.72  ms. 

As  a reasonable  single  correction  when  the  impulse  applied  during 
motion  is  less  than  the  total  impulse,  let  the  impulse  be  that  at 
tj.  = 4.72  ms;  that  is,  I = 1.5  psi-s.  Then 

Wj.  = 1.44  in.  t^.  = 3.95  ms 

For  the  WES  bunker  the  predicted  deflection  is  2.85  inches,  whereas 
the  experimental  deflection  is  0.5  inch;  thus,  the  prediction  is  too 

high.  The  method  therefore  has  to  be  examined  to  account  for  the  high 

deflection  prediction.  The  approximation  that  definitely  leads  to  high 
deflections  is  the  replacement  of  the  blast  pulse  by  a rectangular  pulse 
with  the  same  peak  pressure  and  impulse. 

A less  approximate  method  of  choosing  the  rectangular  pulse  is 
first  to  consider  the  experimental  pulse  shape  as  an  initial  spike 
supposed  on  an  exponential  pulse.  The  impulse  associated  with  the 
pressure  spike  is  only  about  5%  of  the  total  impulse.  Hence,  we  con- 
sider the  whole  pulse  as  an  exponential  pulse  with  the  same  total  im- 
pulse. We  determine  the  peak  pressure  and  decay  parameters  and  a by 
matching  the  experimental  impulses  at  two  different  time  points.  This 
new  peak  pressure  is  used  for  the  rectangular  pulse;  the  impulse  is 
still  the  same.  Mathematically,  the  pulse  may  be  expressed  as 


so  the  corresponding  impulse  is 

! = I.d  ' e'Qt> 

where 
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and  is  the  total  impulse.  For  1^  = 1.8  psi-s  and  I = 1.0  psi-s  at 
time  t = 2 ms,  we  have  a = 102  s”1  and  = 730  psi.  The  above  de- 
flection prediction  procedure  then  gives 

t^  = 2.46  ms,  p^  = 487  psi,  X = 1.41,  = 0.36 

t^  = 3.66  ms 
Wf  = 0.97  in. 

As  a single  correction,  when  the  impulse  applied  during  motion  is  less 
than  the  total  impulse,  let  the  impulse  be  that  at  t^  = 3.66  ms;  that 
is,  I = 1.4  psi-s.  Then  the  central  deflection  and  motion  duration 
are 


= 0.58  in.  t^  = 2.83  ms 

The  above  procedure  for  a more  reasonable  choice  of  the  rectangular 
pulse  results  in  a predicted  deflection  of  1.16  inches  for  the  WES 
bunker,  which  is  much  closer  to  the  experimental  value  of  0.5  inch  than 
the  previous  prediction  of  2.85  inches. 

PRESSURE- IMPULSE  CHARACTERIZATION 

The  analysis  of  Appendix  A treats  the  structural  response  of  a 
clamped  rigid-plastic  beam  to  a rectangular  pressure  pulse  with  a 
symmetric  parabolic  distribution  about  the  center;  the  pressure  is 
zero  at  the  center.  The  interaction  pressure  is  the  spatial  average  of 
the  distributed  peak  pressures.  For  the  1/10-scale  bunker  of  Figure  5, 
the  static  interaction  pressure  is  ”pg  = 344  psi.  The  analysis  is 
applicable  up  to  a pressure  parameter  value  of  X ^ = 5.27,  which  corresponds 
to  an  interaction  pressure  of  p^  = 1810  p9i.  According  to  the  assumption 
of  the  roof  peak  load  being  the  same  aa  the  ground  surface  load,  the 
analysis  is  applicable  up  to  a ground  surface  pressure  of  P = 2715  psi; 
the  corresponding  range  and  impulse  are  11  feet  and  1.1  psi-s  at  1/10-scale 
(Figures  6 and  7).  The  duration  of  the  equivalent  rectangular  pulse 
is  0.4  msec,  so  that  the  loading  is  highly  impulsive  when  X = 


We  rewrite  the  deflection  formula  (14)  as 


W = ‘ [l  + ^ Am  P 2 { 2 - (1  - p)2}] 

2up  X [ 3 s J 


(17) 


where  I = p,t  and  p is  written  instead  of  p . When  X = X , formula  (17) 
d d 1 l 

may  be  written  as 


W 


2 


, 4 , 

+ 3 A2msP2 


2 (2 


(1  - P2)2)] 


(18) 


where  p is  the  location  of  the  plastic  hinge  when  X = X^.  We  choose 
W^  as  the  central  deflection  required  for  structural  failure  to  determine 
I^,  which  now  plays  the  role  of  the  ideal  impulse  Iq  in  the  PI  character- 
ization. The  isodamage  or  structural  PI  curve  is  obtained  by  letting 
W = W^.  Thus,  equating  formulas  (17)  and  (18)  gives 


W2  [ 

1 + ~ x m p2 { 2 - (1  - p )2> 

3 2 s 2 2 

L 

x2P2 

j"l  + | Xmgp2  { 2 - (1  - P)2} 

(19) 


where  p is  a function  of  X. 

Figure  12  shows  the  structural  curve  of  formula  (19)  along  with 
the  universal  structural  PI  curve  of  formula  (1)  where  1^  takes  the 
part  of  the  ideal  impulse  Iq. 


32 


SUMMARY  OF  PROCEDURE 


The  procedure  developed  above  for  obtaining  the  PI  characterization 
of  a shallow-buried  box  bunker  subjected  to  a surface  blast  wave  consists 
of  the  following  steps: 

* Neglect  the  influence  of  the  blast  wave  transit  time. 

* Redistribute  the  pressure  from  a uniform  distribution  on  the 
ground  surface  to  a parabolic  distribution  on  the  bunker  roof, 
having  zero  pressure  at  the  center  and  the  maximum  pressure 

at  the  bunker  edges  to  represent  soil-structure  interaction. 

The  parabola  is  determined  by  assuming  that  the  total  load 
on  the  roof  equals  the  ground  surface  load  over  the  roof. 

(Soil  wave  transmission  and  reflection  effects  are  neglected 
at  this  stage  of  the  development.) 

* Preserve  the  pulse  duration  and  shape. 

* Calculate  the  average  interface  peak  pressure  acting  on  the 
entire  roof  using  the  assumed  parabolic  distribution,  and  set 
it  equal  to  the  surface  peak  pressure;  this  step  equates  the 
total  interface  force  on  the  bunker  and  the  ground  surface 
force  immediately  above,  which  is  consistent  with  static  test 
and  static  finite-element  code  results. 

* Calculate  the  average  peak  interface  pressure  acting  on  the 
roof  slab  (internal  dimension  is  the  slab  span)  and  define 

this  pressure  as  the  interface  pressure  for  PI  characterization. 

* Calculate  similarly  the  average  interface  impulse  on  the  roof 
span  and  define  this  impulse  as  the  interface  impulse  for 

PI  characterization;  this  calculation  is  based  on  the  assumed 
preservation  of  pulse  shape. 

* Neglect  influence  of  soil  wave  transit  times  and  assume  that  the 
lateral  pressure  acts  on  the  entire  sides  of  the  bunker  only 
while  the  roof  pressure  is  acting  (roof  and  side  pressures  in 
phase)  . 

* Calculate  the  lateral  pressure  and  hence  the  roof  thrust  by 
applying  the  Mohr-Coulomb  criterion  with  the  low  cohesion 
value  neglected  and  using  the  surface  pressure  for  the  vertical 
stress;  for  this  step,  the  soil  friction  angle  is  required. 

It  is  assumed  that  the  total  lateral  load  from  the  interaction 
pressure  is  the  same  as  the  lateral  load  away  from  the  bunker 
walls,  where  the  pressure  distribution  is  uniform  (because  of 
neglect  of  wave  transit). 


Derive  the  moment- thrust  yield  condition  for  the  reinforced 
concrete  roof  slab  design;  the  yield  condition  includes  the 
extremes  of  fully  plastic  moment  and  fully  plastic  thrust. 


* Calculate  the  roof  thrust  corresponding  to  the  peak  surface 
pressure,  divide  it  by  the  fully  plastic  thrust,  and  determine 
the  ratio  of  the  resistive  moment  to  the  fully  plastic  moment 
from  the  yield  condition;  hence,  determine  the  resistive  moment. 

* Assume  that  the  roof  is  a clamped  rigid-plastic  beam. 

* Calculate  the  average  static  collapse  pressure  for  the  same 
parabolic  distribution  as  for  dynamic  loadings;  define  this 
pressure  as  the  static  collapse  pressure  required  in  the  PI 
characterization. 

* Replace  loading  pulse  shape  by  a rectangular  pulse  having  the 
same  peak  pressure  and  impulse.  (An  impulse  reduction  is 
possible  if  the  motion  duration  turns  out  to  be  less  than  the 
pulse  duration.) 

* Calculate  the  final  central  deflection  from  the  results  of 
the  dynamic  rigid-plastic  theory  (available  for  rectangular 
pulses  with  pressures  over  five  times  the  static  collapse 
pressure) . 

* Calculate  the  static  collapse  pressure  and  an  approximate  ideal 

impulse  for  P and  I in  the  PI  characterization, 
o o 

* Derive  the  structural  PI  curve  (see  Section  3). 

* Using  PQ  and  IG  plot  the  interface  PI  curve  with  range  as  an 
intrinsic  coordinate. 
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4. 


BURIED  CHARGE  LOADING 


PROCEDURE 

We  have  used  a procedure  similar  to  that  described  in  Section  2 
to  obtain  a Pressure-Impulse  (PI)  characterization  of  the  WES  shallow- 
buried  box  bunkers  subjected  to  loading  from  buried  charges.  Figure  2 
(Section  1)  shows  the  experimental  configuration. 

The  procedure  adopted  includes  the  following  steps: 

Generate  the  interface  pressure  and  impulse  values 
from  formulas  that  fit  the  WES  measurements. 

* Model  the  bunker  panel  next  to  the  buried  charge  as 
a clamped  rigid-plastic  beam. 

* Assume  that  the  pressure  distribution  along  the  beam 
is  parabolic  with  zero  pressure  at  the  center;  define 
the  average  pressure  as  the  interaction  pressure  for 
PI  characterization. 

* Neglect  thrust  because  of  the  proximity  of  the  ground 
free  surface. 

Determine  theoretically  the  static  collapse  pressure, 

P0  , of  a clamped  rigid-plastic  beam  with  the  parabolic 
pressure  distribution. 

* Substitute  the  results  P^,  1^  of  the  WES  incipient  failure 
experiment  and  the  static  collapse  pressure  P0  in  formula  (2) 
to  determine  the  ideal  impulse  I0 

* Plot  the  universal  PI  curve  given  by  formula  (1);  this 
curve  will  pass  through  the  points  P , I . 

* Plot  the  interface  PI  curve  with  range  as  an  intrinsic 
coordinate. 

More  experimental  data  and  at  least  a simple  theoretical  treatment 
for  buried  charges  are  required  to  establish  whether  the  resulting  PI 
characterization  can  be  used  to  predict  failure  loads  for  other  bunker 
designs,  soils,  and  explosive  weights. 
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PI  CHARACTERIZATION 


Free  field  pressures  and  impulses 
may  be  approximated  by  formulas  of  the 


from  buried 
form 


spherical  charges 


I = k.X 

i 


X = R/W 


1/3 


(14) 


with  R as  the  range  and  W as  the  explosive  weight.  Empirical  fits  of 
the  WES  data  are 

”2  1 1 

P = 2650X-4'98  psi  I = 3202X-  ' psi-ms 


(15) 


For  the  incident  pulse  where,  in  X,  the  range  R is  given  in  feet  and 
the  TNT  explosive  weight  W is  given  as  21  lb.  Fits  of  the  interface 
pressure  and  impulse  measurements  are 

P = 2800X-2’45  psi  I = 4444X-  ’ psi-ms 


(16) 


The  parabolic  distribution  of  pressure,  Eq . (5),  with  zero  central 
pressure,  is  taken  to  represent  the  soil-structure  interaction.  The 
fully  plastic  moment  without  thrust  is  Mq;  hence,  the  average  pressure 
required  to  cause  static  collapse  is 

8M 

o 


where  the  beam  span  is  29 . 

In  the  WES  experiment  the  box  bunkers  were  twice  the  size  of  the 
one  shown  in  Figure  5-  Hence,  the  results  in  Section  3 give  the  fully 
plastic  moment  as 

M = 16,400  lb.  in. /in. 
o 

so  that  Eq.  (17)  gives  the  static  collapse  pressure  as 

P = 228  psi 
o 
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In  the  WES  experiments  incipient  structural  failure  occurred  when 
the  range  was  4 feet,  and  catastrophic  failure  occurred  w ien  the  range 
was  2-3/4  feet.  These  results  are  shown  schematically  in  Figure  13. 
According  to  Eq.  (16),  the  interface  pressure  and  impulse  values  at  a 
range  of  4 feet  are 

= 1127  psi  1^  = 2920  psi-ms 

These  values  and  the  static  collapse  pressure  substituted  in  Eq.  (2) 
give 


I = 2329  psi-ms 
o 

Thus  the  universal  structural  PI  curve  can  be  determined  and  plotted 
as  shown  in  Figure  14. 

Elimination  of  X from  Eq.  (16)  and  introduction  of  P and  I give 

o o 

the  interaction  PI  curve 


This  curve  is  also  plotted  in  Figure 

Figure  14  is  the  required  PI  characterization  of  the  box  bunkers. 

As  indicated  above,  further  research  is  necessary  to  establish  whether 

the  characterization  can  be  generalized.  For  example,  the  procedure 

for  determining  the  effects  of  other  explosive  weights  is  simply  to 
1 / 3 

note  that  X = R/W  = 1.45  (R  = 4 feet,  W = 2 lb  TNT)  at  incipient 
1/3 

failure  so  R = XW  gives  the  new  range.  However,  for  larger  charges 
the  interface  PI  curve  may  not  be  the  same  because  of  the  influence  of 
the  free  surface. 
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(a)  VIEW  FROM  CHARGE  <b>  CENTRAL  CROSS 

SECTION 

2-3/4-FOOT  RANGE 


(a)  VIEW  FROM  CHARGE  <b>  CENTRAL  CROSS 

SECTION 


4-FOOT  RANGE 

MA-6203-5 


FIGURE  13  BOX  BUNKER  DAMAGE  FROM  21-POUND  TNT  BURST 
AT  2-3/4-FOOT  AND  4-FOOT  RANGES 
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Appendix  A 


RESPONSE  OF  A CLAMPED  BEAM 
TO  A SYMMETRICALLY  DISTRIBUTED  RECTANGULAR  PULSE 

INTRODUCTION 

The  purpose  of  this  analysis  is  to  determine  the  permanent  contral 
deflection  of  a clamped  beam  caused  by  a rectangular  pulse  that  is 
symmetrically  distributed  about  the  center.  Because  we  are  investigating 
structural  failure,  the  deflections  are  much  larger  than  elastic  deflections. 
Consequently,  the  material  behavior  may  be  considered  rigid-perfect ly  plastic. 
The  choice  of  a rectangular  pulse  greatly  simplifies  the  analy  sis  and  should 
cause  deflections  similar  to  those  caused  by  a blast  pulse  when  the  duration 
and  blast  characteristic  time  are  a fraction  of  the  structural  response 
time.  The  symmetrical  pressure  distribution  is  an  approximate  representation 
of  soil-structure  interaction  and  is  specialized  to  a parabolic  distribution 
having  the  minimum  pressure  at  the  center. 

The  beam  supports  are  allowed  to  approach  each  other  so  that  negligible 
tensile  membrane  forces  develop.  Also,  while  the  pressure  pulse  is  acting 
it  is  assumed  that  a constant  thrust  is  applied,  having  as  its  only  effect 
the  creation  of  an  increased  fully  plastic  resistive  moment  in  the  beam. 

When  the  pressure  pulse  ends,  the  thrust  also  ends. 

In  rigid-plastic  theory,  the  pulse  pressure  must  exceed  the  static 
collapse  pressure  to  produce  permanent  deformation.  The  static  and  dynamic 
pressures  therefore  differ  in  magnitude,  but  the  distributions  are  chosen 
to  be  the  same. 

STATIC  COLLAPSE  PRESSURE 

As  the  static  pressure  p(x)  is  increased,  plastic  hinges  appear  at 
the  center  and  at  the  supports  to  form  the  static  collapse  mechanism 
(Figure  Al).  At  each  hinge,  the  fully  plastic  moment  is  M^  and  the  equa- 
tion of  equilibrium  at  incipient  collapse  is 


A- I 


1 


i 


J 


(b)  THREE-HINGED  STATIC  COLLAPSE  MODE 


» 1 

R 


(c)  FORCES  AND  MOMENTS 

MA-6203-1 

FIGURE  A1  STATIC  COLLAPSE  OF  A RIGID-PLASTIC  CLAMPED  BEAM 


(Al) 


where 


F (i)  = 2 M 
s s 


l 

I 


FsU)  = | (l  - x)  ps(x)  dx 


( A2) 


L 


A parabolic  distribution  of  pressure  is  described  by 

2 


p ■ + q 


2 „2 


so  the  average  pressure  pg  and  Fg(^)  from  (A2)  are 


ps  ' % + 3 "2 


- ? (%  + £ ’2)  l2 


and  the  equilibrium  equation  (Al)  becomes 

, 4 M 

. 1 s 

qo  6 q2  = 


(A3) 


(A4) 

(A5) 


(A6) 


For  zero  central  pressure  we  have  q^  = 0;  hence,  Eqs.  (A4)  and  (A6) 
give  the  static  collapse  pressure 

8 M 


Ps  = 'Yr 


(A7) 


DYNAMIC  RESPONSE  IN  MECHANISM  1 

If  the  suddenly  applied  pressure  is  only  slightly  greater  than  the 
static  collapse  pressure,  it  is  reasonable  to  postulate  that  the  beam 
will  deform  in  the  static  collapse  mode  because  the  inertial  forces  are 
not  large  enough  to  affect  the  mechanism.  Consequently,  the  displacement 
field  is  taken  as 


where  W(t)  is  the  central  deflection. 


A- 3 


In  the  notation  of  Figure 


A2,  the  equations  of  motion  are 


22 

dx 


-p  + M> 


(A9 


when  rotary  inertia  is  neglected. 

We  analyze  the  motion  in  two  phases.  Phase  1 is  the  initial  loaded 
phase  and  Phase  2 is  the  subsequent  unloaded  phase. 


Phase  1 p = pd(x),  t < td 

When  we  substitute  the  displacement  field  (A8)  in  Eqs . (A9),  integrate 
spatially,  and  use  the  hinge  conditions  M(o)  = Md  and  M(.£)  = -Md>  we 
obtain  the  governing  equation 


where 


W 


Fd(*> 


Fd(£)  = j (£  - x)  pd(x)  dx 


(A10) 


(All) 


Dots  over  symbols  indicate  time  differentiation. 


We  introduce  a pressure  parameter  defined  by 

PjW 

- ps(x> 

so  that  from  Eqs.  (Al),  (A2) , and  (All)  we  have 


( A12) 


Fd(f)  = KFs(£)  = 2 MsA 
Eq.(AlO)  may  therefore  be  written  in  the  form 


W 


( Am 

s 


m ,) 
d 


( Al  3) 


(A14) 
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(a)  DYNAMICALLY  LOADED  CLAMPED  BEAM 


w = W(t)  (I-*) 


(b)  MECHANISM  1 


M ( M + \ 


Q f q + ?7> 



'^2Wdx 

(c)  FORCES  AND  MOMENTS  ON  RIGID  BEAM  ELEMENT 

MA-6203-3 

FIGURE  A2  DYNAMIC  COLLAPSE  OF  A RIGID-PLASTIC  CLAMPED  BEAM 


where 


M 

s 

m = TT~ 
s M 

o 


!*d 

md  M 


( A15) 


The  beam  central  velocity  and  deflection  at  time  t.  are  therefore 

d 


6 M 


Wd  " 


7T  (Am 
2 s 


3 M 


mdHd 


Wd  = 


M * 


y (Am  - m )t3 

2 s d d 


(A16) 


Phase  2 


p = o,  t > t 


We  assume  that  in  this  phase  the  thrust  is  no  longer  acting! 

so  that  the  fully  plastic  moment  decreases  instantaneously  at  time  t, 

d 

from  Md  to  Mq.  Setting  pd>  and  hence  A,  to  zero  in  Eq . (A14)  and  replacing 

M,  by  M give 
d o 


6 M 


W = - 


( A17) 


VSL 


Temporal  integration  of  Eq . (A17) , substitution  of  the  initial  velocity 
condition  ( A16) , and  setting  W(f^)  = 0 lead  to  the  equation  for  the  dura- 
tion of  motion  tf. 


1 + Am  - m. 
s d 


(A18) 


Another  temporal  integration  and  substitution  of  the  initial  dis- 
placement condition  (A16)  gives  the  final  central  deflection  as 

\ 

2 

3 M td 

W,  = — " — l Am  - m,)(l  + Am  - m, ) 
t .2  s d s d 

P*. 


(A19) 
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r 


l 


For  a parabolic  pressure  distribution  with  zero  central  pressure 
we  have  the  average  pressure 


8 M 


Pd  = ^Ps 


Xm 


( A20) 


so  that  formula  (A19)  for  the  final  central  deflection  becomes 


W 


f 


3 Vd2 


8p 


(1  + \m  - m , ) 
s d 


(A21) 


Pressure  Limit 

Pressures  considerably  higher  than  the  static  collapse  pressure 

initiate  mechanisms  different  from  mechanism  1 because  of  the  greater 

importance  of  the  inertial  forces.  The  upper  limit  of  p.  or  X to  initiate 

d 

mechanism  1 is  determined  by  a change  of  sign  at  a central  or  support 
plastic  hinge  of  the  second  spatial  derivative  of  the  moment  because  this 
implies  incipient  violation  of  the  yield  condition.  From  Eqs.  (A9) 


9 

a v 

pd  + 11  a 2 

t 


( A22) 


After  substituting  the  displacement  field  ( A8)  in  ( A22) and  setting  in  turn 
the  hinge  locations  x = o and  x = 1,  we  see  that  the  upper  limit  of  X is 
obtained  by  a sign  change  at  the  center.  This  limit  for  zero  central 
pressure  is  X^  = 1;  hence,  mechanism  1 cannot  be  initiated  by  this  special 
parabolic  distribution  of  pressure. 

DYNAMIC  RESPONSE  IN  MECHANISM  2 

For  load  parameters  slightly  higher  than  X^  the  beam  responds  in 
mechanism  2 (Figure  A3), which  has  four  plastic  hinges.  While  the  rectan- 
gular pulse  is  acting, the  hinges  remain  stationary  at  x = x^  and  the 
central  portion  of  the  beam  moves  by  rigid  body  translation. 
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Phase  1 


The  unknowns  in  Eqs . (A24) , (A25) , and  (A28)  are  the  central  deflec- 
tion W,  the  central  moment  M(o),  and  the  hinge  location  x^.  The  hinge 
location  can  be  determined  from  Eq.  (A28) . For  the  parabolic  pressure 
dis tribution  ( A3)  with  q^  = 0,  Eq . (A28)  reduces  to 


From  Eq.  (A24)  the  central  velocity  and  deflection  at  time  t^  are 

- 


PdCd 


Wd  = tx 


W 

d 2 M- 


(A30) 


where,  for  zero  central  pressure 


Pd  = Xps 


8 M 


T^l 


(A31) 


Phase  2 


p = 0,  ^ < t < t 


After  the  pressure  has  been  applied,  the  hinges  travel  toward  the 
center  of  the  beam;  let  the  arrival  time  be  t^.  The  velocity  field 
corresponding  to  this  phase  of  deformation  is 


<3w  _ j 

at  ' ^ 


i - X 


i 


0 s x £ x (t) 


x^(t)  £ x £ SL 


( A32) 


where  V = V(t)  is  the  velocity  of  the  center  of  the  beam. 

In  the  region  0<  x<  x^(t),  the  velocity  field  (A32) , the  equations 
of  motion  ( A9 ) , and  initial  conditions  (A30), give  at  time  t 


V = w = w , 
c c d 


W = W,  + W,  (t  - t ,) 
c d d d 


( A33) 


In  the  region  x^t)  < x < l,  the  velocity  field  (A32),  the  equations  of 
motion  (A9),  initial  conditions  (A30),and  moment  conditions  MCx^)  = Mq 
and  M(J.)  = -M  , give  at  time  t 


t - t , 
c d 


(A34) 
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r 

Phase  3 P = 0,  < t< 

The  remaining  deformation  occurs  in  mechanism  1,  which 

collapse  mode,  arid  has  the  velocity  field 

is  the  static 

(A35) 

The  initial  conditions  of  this  phase  are  Eqs . (A33) . The  solution  of  the 
equations  of  motion  (A9)with  velocity  field  (A35) and  initial  conditions 

is 

wf  -"c-K (h  - 0 

( A36) 

M.V  i2 
c 

tf  ' tc  6 M 

0 

( A37) 

where  is  the  required  final  central  deflection  and  t^  is 

of  motion.  Explicitly,  the  final  central  deflection  is 

the  duration 

( A38) 

For  the  parabolic  pressure  distribution  (A3) with  = 0,  the  initial 

hinge  location  is  the  solution  of  Eq.  (A29),and  the  central  deflection 

is 

Wf  = 2 n P*  | 1 + 3 VT  [2  ‘ (X  - Pl)  J j 

(A39) 

and  the  duration  of  motion  is  given  by 

^ - 1 + § vi  [3  - (‘  -pi)2_ 

(A40) 

A- 11  j 

1 

/ 

L — inn  ; — J 


r 


"1 


Pressure  Limit 

Mechanism  2 requires  that  the  central  bending  moment  lie  in  the 
range  -M^ < M(o)  < . The  largest  value  of  X that  allows  this  condition 

to  be  satisfied  is  determined  by  Eq . (A25)  with  M(o)  = -M^ , that  is. 


(A41) 


For  a parabolic  pressure  distribution  with  zero  central  pressure,  formula 
(A41)  gives 


A „m 
2 s 


_4 


(A42) 


and  Eq. (A29)  determines  by  becoming 


Pi  + 5 


0 


(A43) 


Thus,  to  find  \2>  we  solve  Eq.  (A43)  for  and  substitute  in  Eq . (A42) . 
Calculations  show  that  = 0.72  and  A^  = 5.27.  The  pressure  limit  for 
activating  mechanism  2 is  therefore  over  five  times  the  static  collapse 
pressure. 


■ 
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Appendix  B 

REINFORCED  CONCRETE  BEAM  YIELD  CONDITION 


INTRODUCTION 

This  appendix  outlines  the  derivation  of  the  yield  condition  for 
the  reinforced  concrete  slabs  of  the  box  bunker.  The  yield  condition 
is  a relationship  between  the  bending  moment  M and  the  axial  thrust  N 
when  the  cross  section  yields  plastically  over  its  entire  surface.  In 
the  derivation  the  following  assumptions  are  required: 

* Plane  sections  remain  plane  during  deformation 

* Perfect  bonding  exists  between  the  steel  and  concrete 

* The  steel  is  perfectly  plastic  (no  hardening) 

* The  concrete  is  perfectly  plastic  in  compression  and  has 
no  tensile  strength 

* Round  reinforcing  bars  may  be  approximated  by  square  ones 

MOMENT  AND  THRUST  FORMULAS 

Figure  B1  shows  the  geometrical  nomenclature  of  the  slab;  the 

usual  round  bars  have  been  replaced  by  square  ones  having  the  same 

area  for  algebraical  convenience  only.  Figure  B2  shows  the  four  stress 

distributions  that  can  arise  for  the  bunker  slab  cross  section.  Figure  B2a 

applies  when  the  neutral  axis  (plane  of  zero  strain)  at  depth  c intersects 

the  upper  steel,  that  is,  when  d'  < c < d'.  For  the  properties  of  the 

bunker  slab,  the  limiting  conditions  of  zero  thrust  and  fully  plastic 

moment  M fall  in  this  case.  As  the  thrust  is  increased  and  the 
o 

moment  is  changed  to  keep  the  section  fully  plastic,  the  depth  of  the 
neutral  axis  increases,  resulting  in  the  stress  distributions  of 
Figures  B2b,  c,  and  d.  The  limiting  condition  is  the  fully  plastic 


B-l 


thrust  N with  no  resistive  moment.  Formulas  for  the  moments  and  thrusts 
o 

are  derived  below  for  the  four  cases  in  Figure  B2.  Dimensionless 
quantities  that  are  convenient  for  expressing  results  are 


a 


(Bl) 


Case  1:  d^  < t < d^ 

The  stress  distribution  in  Figure  B2a  may  be  represented  by  the 
following  forces  and  corresponding  moments  (moments  about  section  mid- 
height) : 


The  compressive  axial  force  and  the  binding  moment  are 

-N'  + N - N ' - N 

c su  si  s 

- M ' + M ' - m'+M 

C SU  S 1 S 

(s  - 1)6  ly  - 6 j \ - 36^2  - vj-  set  (B2) 


N = N ' 
c 


M = M 


Hence 


bho 


= y + 


B-4 


which  gives  the  position  of  the  neutral  axis  for  bending  alone  provided 

6'  < 7q  < 6^.  For  the  bunker  design  of  Figure  B2a  this  inequality  is 

satisfied.  Thus,  setting  y = y in  Eq.  (B3)  gives  the  fully  plastic  moment 

o 

M without  thrust, 
o 


Case  2 : d ' < c < d 

2 1 

The  stress  distribution  in  Figure  B2b  may  be  represented  by  the 
following  forces  and  corresponding  moments: 


The  compressive  axial  force  and  bending  moment  are 

N = N ' - N ' + N ' - N 
cess 

M = M ' - M ' + M ' + M 
c c s s 


2M 


2 

bh  p 


= y (1  - V>  + (2s  - 1)  fy  (1  - 26') 


(B6) 


Case  3rd  < c < d 
1 2 

The  stress  distribution  in  Figure  B2c  may  be  represented  by  the 
following  forces  and  corresponding  moments: 


N'  = be 


N'  = a a 
c u 


N'  = a a 

S y 

N = a (c  - d )CX 
su  1 y 


N = a(d  - c )cr 
si  2 y 


N = a(c  - d )(7 

cu  1 u 


M 


M 


- = N'/i  - h) 

c c \2  2/ 

' =N'(*-d') 

B c \2  / 

= N [l(c  + d ) - -1 
su  su  L 1 2j 

= N i(c  + d ) - - 1 

3 1 sl  L 2 2 J 

= N [i(c  + d ) - - 1 
cu  cu  L 1 2 J 


The  compressive  axial  force  and  bending  moment  are 

n = n'-n,  + n'  + n -n  -n 
c c s su  sl  cu 


M = M'  - M'  + M'  - M +M  +M 

c c s su  sl  cu 


Hence , 


= y + (s  - l)  @ (y  - 61>  - S3  (62~  y)  + (s  - l)* 


(B7) 


2M 


2 

bh  o 


= y (1  - y)  - (s-l)3  (y-61>  £(y  - 6i>  - lj  + Sg  (62~  y) 


+ ( s -1  )fy  (1 
(B8) 
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V’1] 

- 26') 


Case  4: 


d < c < h 
2 


The  stress  distribution  in  Figure  B2d  may  be  represented  by  the 
following  forces  and  corresponding  moments: 


The  compressive  axial  force  and  bending  moment  are 


N = N'  - N'  + N'  + N - N 
c c s s c 


M=M'-M  +M'-M  + M 

c c s s c 


Hence, 


N 

= y + 2 (s  - l)a 

(B9) 

bhCT 

u 

2M 

= y (i  - y) 

(BIO) 

O 

bh  CT 

u 


When  y = 1,  the  section  sustains  only  an  axial  force  N , which 

o 

from  Eq,  (B9)  is 


= 1 + 2 (s  - l)a  (Bll) 

bhCT 

u 
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Dimensionless  forms  of  the  moment  and  thrust  that  are  convenient  for 


expressing  the  yield  condition  are 

m = M/M  n = N/N  (B12) 

o o 

where  Mq  and  NQ  are  given  by  Eq . (B3)  with  y = y and  Eq . (Bll).  The  simplest 
way  to  construct  the  yield  curve  is  to  regard  the  dimensionless  depth 
of  the  neutral  axis  as  a parameter  and  to  call  upon  the  moment  and 
thrust  formulas  according  to  the  range  in  which  y falls. 


SLAB  YIELD  CONDITION 

Data  corresponding  to  the  design  in  Figure  5 are: 


Slab  thickness 
Spacing  of  reinforcing  bars 
Diameter  of  reinforcing  bars 
Area  of  reinforcing  bar 
Side  of  equivalent  square  bar 
Depth  to  center  of  upper  steel 
Depth  to  center  of  lower  steel 
Depth  to  top  of  upper  square  bar 
Depth  to  bottom  of  upper  square  bar 
Depth  to  top  of  lower  square  bar 
Depth  to  bottom  of  lower  square  bar 
Crush  strength  of  concrete 
Yield  strength  of  steel 


h = 2.8  in 


b = 2.0  in 


D = 1/4  in 

2 

A = 0.049  in 
s 

a = 0.222  in 


d ' = 0.4  in 


d = 2.4  in 

d'  = 0.289  in 

d ' = 0.511  in 
2 

dj^  = 2.289  in 

d2  = 2.511  in 

O = 6,000  psi 
u 

CT  = 72,000  psi 

y 


From  these  data  the  following  parametric  values  were  obtained: 

2 

C*  = a /bh  = 0.00877 
0 = a/b  = 0.111 

g ' = d7h  = 0.103 

6 ' = d/h  = 0.143 

6 ' = d'/h  = 0.182 
2 2 

6i  = d^/h  = 0.818 

6 = d/h  = 0.857 

6 = d /h  = 0.897 

2 2 

s = 0/0  =12 

y u 

The  resulting  fully  plastic  moment  with  no  thrust,  fully  plastic 
thrust  with  no  moment,  and  the  dimensionless  depth  of  the  neutral  axis 
when  the  moment  acts  alone  are 

M = 4,100  lb.  in/in 
o 

N = 20,000'  lb/in 
o 

y = o.i34  < y < 

o \ 1 o 2/ 

Table  B1  lists  the  values  of  m and  n generated  by  means  of  the 

parameter  y . Figure  B3  is  the  resulting  yield  curve.  An  interesting 

feature  of  the  curve  is  that  when  the  thrust  is  about  41%  of  the  fully 

plastic  thrust  N , the  resistive  moment  is  2.26  times  the  fully  plastic 

moment  M . 

o 
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Table  B1 


YIELD  CURVE  VALUES 


Y Range  and 

Formulas 

Y 

m 

n 

6;  < y < s 2 

0.134 

1.00 

0 

0.14 

1.10 

0.019 

(B2)  (B3 ) 

0.16 

1.38 

0.079 

0.18 

1.65 

0.138 

^2  K y 61 

0.20 

1.74 

0.160 

0.30 

2.03 

0.244 

0.40 

2.20 

0.328 

(B5)  (B6) 

0.50 

2.26 

0.412 

0.60 

2.20 

0.496 

0.70 

2.03 

0.579 

0.80 

1.74 

0.663 

61  < V 62 

0.82 

1.65 

0.685 

0.84 

1.38 

0>745 

(B7)  (B8) 

0.86 

1.10 

0.804 

0.88 

0.80 

0.863 

62  < Y < 1 

0.90 

0.516 

0.916 

0.94 

0.323 

0.950 

(B9)  (BIO) 

0.98 

0.112 

0.983 

1.00 

0 

1.0 
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